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Abstract 

The lower-dimensional maximal polytopes associated with an iteration stable 
(STIT) tessellation in M"^ are considered. They arise in the spatio-temporal construc- 
tion process of such a tessellation as intersections of {d — 1) -dimensional maximal 
polytopes. A precise description of the joint distribution of their birth-times is ob- 
tained. This in turn is used to determine the probabilities that the typical or the 
length- weighted typical maximal segment of the tessellation contains a fixed number 
of internal vertices. 
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1 Introduction and results 

Random tessellation theory is an active field of current mathematical research. Besides 
theoretical developments, there are numerous applications of random tessellations for 
example in the study of random structures in biology, geology and other sciences; cf. [U 
El [ini [HI US] • Apart from the classical models, such as the well known Poisson hyperplane 
or Poisson- Voronoi tessellations for which we refer to [3l [101 [13 [IS] , the class of iteration 
stable (or STIT) tessellations has attracted considerable interest in recent times; see 
[H [71 [Hi [HI [151 [13 [IB] and the references cited therein. In particular and as discussed 
in [8], the STIT tessellations may serve as a reference model for hierarchical spatial cell- 
splitting and crack formation processes in natural sciences and technology, for example 
to describe geological or material phenomena or aging processes of surfaces. 

Intuitively and within a compact convex polytope W C M"^ (we assume that d > 2 
in this paper) with positive volume, their construction can be described as follows. At 
first, the window W is equipped with a random lifetime. When the lifetime of W runs 
out, we choose a hyperplane H, which divides W into two non-empty sub-polytopes 
and W~ . Now, the construction continues independently and recursively in and W~ 
until some fixed time threshold f > is reached; see Figure [Tj The outcome Yw{t) of this 
algorithm is a random subdivision (tessellation) of W into polytopes (called cells in the 
sequel). To be more precise, we have to specify the lifetime distribution of the cells and 
the law of the cell-separating hyperplanes. For this purpose, let us write Ti for the space 
of affine hyperplanes and Hq for the subspace of linear hyperplanes in M*^. Furthermore, 
let A be a measure on Ti, which admits the decomposition 

(1) f f{H)A{dH)= f [ fiHo + x)e„Adx)QidHo), 
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Figure 1: States of the random cell division process at different time instants; the re- 
spective new segments are dashed. Here, W is a square. 



where / is a non-negative measurable function on Ti, is the Lebesgue measure on 
Hq and Q is a probability measure on T-Lq. We require that Q is non-degenerate in the 
sense that span{n(i?) : H G suppQ} = M'^, where n{H) is the unit normal of H lying 
(for definiteness) in the upper unit half-sphere. For a polytope c write (c) = {H € % : 
cD H ^ 9} for the collection of hyperplanes hitting c. Now, the lifetime of a cell c in 
the above construction is chosen to be exponentially distributed with mean A((c))^^ (if 
Q is the uniform distribution on Hq, A((c)) is just a dimension dependent multiple of 
the mean width of c). Moreover, we choose the hyperplane splitting c according to the 
(conditional) law A( • |(c)). It is exactly this choice, which makes mathematical analysis 
of the so-constructed tessellations possible and completes the description of the above 
algorithm; see [9] for more details. In particular, we notice that the exponential lifetimes 
of the cells ensure that the construction enjoys the Markov property in the continuous 
time parameter t. 

Besides looking at the local tessellation Yw{t) within W, it is convenient to extend 
Yyir(t) to a whole space random tessellation Y{t) in such a way that for any W as above, 
Y{t) restricted to W has the same distribution as the previously constructed Yw{t) (this 
is possible by consistency according to the main result in [9j). We call Y{t) a STIT 
tessellation of since Y{t) enjoys a stochastic stability under iterations as explained 
later (that Y{t) is in fact a tessellation is ensured by the non-degeneracy assumption on 
Q). It is important that the form ([T]) of A implies that the distribution of Y{t) is invariant 
under spatial translations, i.e., the shifted tessellation Y{t) +x has the same distribution 
as Y{t) for any x G M'^. 

With Y{t), a number of geometric objects are associated. To introduce them, we 
write M'P'^^'i family of cell-splitting hyperplane pieces that are introduced during 

the recursion steps in the above algorithm until time t (these are the dashed segments 
in Figure More generally, for k = 0, ... ,d — 2 we denote by MV^^ the family of 

A:-dimensional faces of members of AdV^j^}^^. For k = 0, . . . ,d — 1 we call A^7-'[,*^ the class 
of fc-dimensional maximal polytopes of y(t). They are the natural building blocks of 
Y{t) and its lower-dimensional face-skeletons. We also consider /c-dimensional weighted 
maximal polytopes, where the intrinsic volumes Vj, < j < k, constitute the weights. 
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To define them, fix A; E {0, . . . , d — 1} and j £ {0, . . . ,k}, let m(p) be the circumcenter of 
a polytope p, write for the (measurable) space of fc-polytopes with circumcenter at 
the origin o and let be the centered ball in M*^ with radius r > 0. Now, we introduce 
a probability measure P^*^- on as follows: 

E 1{P C B^} l{p- m{p) E A} Vj {p) 



2) Pg(^)-lim^^ 



E E l{pCB!i}Vjip) 



where A is a Borel subset of (following the proof of Equations (10) and (11) in [12], 
it can be shown that the limit is well-defined; another argument can be given by using 
([9]) below). A random polytope with distribution pg. is called a Vj-weighted typical 

fc-dimensional maximal polytope of Y{t) and will henceforth be denoted by MPg-. 
If j = 0, this is the typical A;- dimensional maximal polytope and for j = k we obtain 
the volume-weighted typical fe-dimensional maximal polytope of the STIT tessellation 
Y{t), which are two classical objects considered in stochastic geometry; cf. [31 [T2l [T3]. 
For example, MP^*q is the typical maximal segment, whereas MP^*]^ is the length- 
weighted typical maximal- segment. However, our approach here is more general 
and interpolates between the typical {j = 0) and the volume- weighted typical (j = k) 
maximal polytope of dimension k. 

Any A;-dimensional maximal polytope of Y{t) is by definition the intersection of d — k 
maximal polytopes of dimension d — 1. In view of the spatio-temporal construction 
described above, each of these polytopes has a well-defined random birth-time. We denote 
these random variables by /3i, . . . , f3d-k and order them in such a way that < /3i < . . . < 
Pd-k < t holds almost surely. Our first result describes the joint distribution of these 
birth-times; in the special cases d = 2ord = 3, k = l and j = 01 j = 1, the formula is 
known from (TUZKIH]. 

Theorem 1. Given d > 1, E {0, . . . ,d — 1} and j E {0, . . . ,k}. The joint distribution 
of the birth-times /3i, . . . , Pd-k of the Vj -weighted typical k-dimensional maximal polytope 
of the STIT tessellation Y{t) has density 

k—j 

(si, . . . , Sd-k) ^ {d- j){d-k- 1)!^ 1{0 < si < . . . < Sd-k < t} 

with respect to the Lebesgue measure on the {d — k)- dimensional simplex {0 < si < . . . < 
Sd-k < i}, which is independent of the hyperplane measure A. In particular, if j = k we 
obtain the uniform distribution. 

We turn now to an application of Theorem [l| where we consider the typical and the 
length- weighted typical maximal segment MP^ q and MP^*^]^ of Y{t), respectively. These 
segments may have internal vertices, which arise at the time of birth of the segment 
(when d > 3) and thereafter subject to further subdivision of adjacent cells; see Figure [l] 
for an illustration in the planar case. With the help of Theorem [T] we can determine the 
probabilities pi^o(^) and pi^i(n) that the typical or the length- weighted typical maximal 
segment of Y{t) contains exactly n E {0,1,2,...} internal vertices (we suppress the 
dependency on t in the notation of these probabilities since they are independent of the 
time parameter t). 

Theorem 2. The probabilities pi^o(^) ^.i^'d pi,i(ra) are given by 
Pi, o{n) = d{d - 2)1 / .../ —^7— dsi . . . dsd_i 

Jo Jo Jo [d-t- Sd-l - Sd-2 - ... - Sl)"+^ 
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and 

Pi,i(n) = (n + l)(d-l)! 



JO 



I-S2 s2_^ . ^ _ 2sa_, - Sd-2 - ... - SiY 
td-1 (^d-t- Sd-1 - Sd-2 - ... - Sl)"+2 

dsi . ..dsd-i- 



Moreover, pifi{n) and pi^i(n) are independent of t and A. In the mean, the typical 
maximal segment has ^{d"^ —d+2) / {d—1) internal vertices in dimension d >2, whereas the 
length-weighted typical maximal segment in space dimension d> 3 has (d^ — 2d+4)/((i— 2) 
(the mean is infinite if d = 2). 

In the planar case d = 2, pi^Q{n) is known from [71 [T7], whereas for d = 3 the 
formula for pi o(ra) has been established in |18j by different methods. Our approach in 
the present paper is more general and allows to deduce the corresponding formula for the 
length- weighted maximal segment as well as to deal with arbitrary space dimensions. To 
provide a concrete example, take d = 3 and consider the length-weighted typical maximal 
segment. Here, we have 

63 

pii(O) = 5 + 18 In 2 -— In 3 0.173506, 

657 

pi,i(l) = 28 + 901n2- —In 3^ 0.159712, etc. 

' 8 

The mean number of internal vertices equals 7 in this case. 

The rest of this note is structured as follows. In Section [2] we recall some important 
facts about STIT and also about Poisson hyperplane tessellations, which prepare for the 
proof of our results. These are presented in the two final sections [3] and [4} 



2 Some preliminaries 

Iteration of tessellations. Below we will exploit the fact that the tessellations Y{t) are 
stable under iterations. To explain what this means, let < s < t and define the itera- 
tion Y(s)SY{t) of y(s) and Y(t) as the tessellation that arises by locally superimposing 
within the cells of Y{s) independent copies of Y{t). Formally, let {Yc{t) : c £ Y{s)} be 
a family of i.i.d. copies of Y{t), which is indexed by the cells of Y{s), and which is also 
independent of y(s). Then, 

y(s) ffly(t) := {cn c : c e Y{s), c e Y^t), intcnintc 0}. 

The STIT tessellations are stable under iteration in that the distributional equality 

(3) Y{s)aY{t) = Y{s + t) 

holds for any s,t > 0. In other words, the results are the same in distribution when we 
either run the above cell-division algorithm from time s to time s + t or perform at time 
s an iteration of Y(s) and Y{t); cf. [6l [9]. This will play an important role in the proofs 
below. 



STIT scaling. We collect here two implications of the scaling property of a STIT 
tessellation Y{t). Globally, it says that the dilation tY{t) of Y{t) by factor t has the same 
distribution as Y[l), the STIT tessellation with time parameter 1, i.e., 

(4) tY{t)^Y{l). 
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For a polytope W we also have the local scaling tYw{t) = Ytw{^)', see [9j for example. 
Let us denote by the density of the j-th intrinsic volume of MV^I^ , that is, 



(5) ^(*^.^=lim^E nP^Bf}V,{p), 

where is the volume of Bf (one can in fact show that this limit is well-defined, see 
[H]). Using Q, the definitions ^ of MP^*^- and ([s]) of q'^^- as well as the homogeneity 
of the intrinsic volumes, one easily shows the following two facts. 

Lemma 3. For t > 0, k £ {0, . . . ,d — 1} and {0, . . . ,k} it holds that 



h) ^Vi{M¥f.) = t-^'EVi{MV^^]). 



Besides these two scaling relations, also the exact values of pf^. and EFifMPl*^.) are 

I — I 

known from [M] or can be determined with the help of Lemma below, but they are not 
important for our purposes. 

STIT intersections. Let -E be a /c-dimensional affine subspace of M'^. Then Y[t)r\E is 
a STIT tessellation within E, i.e., the sectional tessellation Y{t) n is also stable under 
iterations. It has the property that the (A; — l)-volume density of its cell boundaries is 
proportional to t with a proportionality constant depending on d, k and the hyperplane 
measure A (this follows from standard intersection formulae for surface processes 
Theorem 4.5.3]). In particular, the intersection of the STIT tessellation Y(t) with a line 
L = spanu, where u £ (upper unit half-sphere), is a Poisson point process with 

intensity tA{{u)) (here u has to be interpreted as the line segment connecting the origin 
with u); cf. |9j. 

Poisson hyperplanes. Let us denote by PHT(t) a Poisson hyperplane tessellation 
in M'' with intensity measure tA; cf. [31 [13] for definition. Now, similarly as for the STIT 
tessellations, for /c G {0, . . . , d — 1} and j € {0, . . . ,k} we denote by ^. the Vj- 
weighted typical A;- face of PHT(t). Formally, its distribution is defined by replacing 

the class MV^^ of A;-dimensional maximal polytopes of Y{t) in ^ by the class of A;- faces 
of the Poisson hyperplane tessellation PHT(t). 

We make use of the scaling property of Poisson hyperplane tessellation, which is 
similar to the scaling property (|4]) of a STIT tessellation. Formally, it says that the 
dilated Poisson hyperplane tessellation tPHT(t) has the same distribution as PIIT(l), 
i.e., 

(6) tPHT(t) = PHT(l). 

This follows directly from the uniqueness theorem for Poisson processes [13', Theorem 
3.2.1] and the form of the intensity measure A; recall ([T]). 

A distributional identity. In our arguments below, we need an identity describing 
the distribution of MP^*^. in terms of weighted faces in Poisson hyperplane tessellations. 
In |141 Theorem 3] this fundamental connection between the STIT and the Poisson hy- 
perplane tessellations has been established for j = by martingale techniques and the 
theory of piecewise deterministic Markov processes. For our purposes we need a slight 
generalization of this identity for arbitrary j G {0, . . . , A;}, the proof of which resembles 
the argument of Lemma 4 in [18| - designed for d = 3, A; = 1 and j = 1. 
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Lemma 4. Let d > 1, A; G {0, . . . , (i — 1}, j G {0, . . . , /c} and t > 0. It holds that 
mMPi>)= /' '"-y E/(Fr7'-')d» 

J 

/or any non-negative measurable function f : — )• M. 

It is worth pointing out that the density s i— )• {d—j)s'^^^^^t^^'^~^^ is the marginal dens- 
ity of the last birth time Pd-k of the Ij-weighted typical maximal polytope of dimension 
k in Theorem [TJ 

3 Proof of Theorem [T] 

Throughout this section we let A be a fixed hyperplane measure with representation ([T]) 
and Y{t) be a STIT tessellation constructed with the hyperplane measure A until time 
t > 0. The plan of the proof of Theorem [l] is first to show the formula for j = k and then 
to establish the formula in full generality from the special case. 

3.1 The case j = k 

Lemma 5. Let d > 1 and k G {0, . . . ,d — 1}. The joint distribution of the birth-times 
Pi, . . . , Pd-k of the volume-weighted typical k-dimensional maximal polytope MP[,*^^ of 
Y[t) is the uniform distribution on {0 < si < . . . < Sd-k < t}; which has density 

{d-k)\^,^ 

(si, . . . , Sd-k) ^ ^d-k 1|0 < si < • • • < Sd-k < t}- 

Proof. Recall that MP^*^- is the intersection oi d — k maximal polytopes of dimension 
d — 1, lying on hyperplanes Hi, ... , Hd-k and having birth-times Pi, ... , Pd-k- We are 
now going to calculate the probability 

(7) P(/3i G {Q,si),...,Pd-k G {sd-k-i^Sd-k)), 

where si < . . . < Sd-k are fixed. In other words, we want to calculate the probability 
that hyperplane Hi is born during the time interval (0, si), hyperplane H2 is born during 
the time interval {si, S2) etc. until hyperplane Hd-k and the time interval {sd-k-i-, Sd-k)- 
To evaluate this probability, we formally mark every member of MV^^ with its asso- 
ciated d — k birth-times Pi, ... , Pd-k- This yields a marked point process on the product 
space c^fc X (0,00)'^"'^, where ^k is the space of /c-dimensional polytopes in M'^. In this 
context, the probability in ([T]) is just the mark distribution of this point process evalu- 
ated at (0,si) X ... X {sd-k-i, Sd-k)- According to the general theory of marked point 
processes (cf. |13i Chapter 3]), this equals ''^''"'"^V^fci) where g^^jj-'-''^'^-'"^^ [g the 

A:-volume density of those fe-dimensional maximal polytopes of Y{t) whose birth-times 
satisfy the constraints Pi G {0, si), . . . , Pd-k G {sd-k-i, Sd-k)- Since = t'^~''g^^j, ac- 
cording to Lemmajsja), it remains to determine £>^'^^' This is done recursively 
using that 

Y{t) = Y{si) ffl Y{s2 -si)m...m Y{sd-k - Sd-k-i) ffl Y{t - Sd-k), 

which is valid by iteration stability of Y{t); recall ([s]). Thus, such a A;-dimensional max- 
imal polytope must be contained in the intersection of d—k processes of (d— l)-dimensional 
polytopes having (d— l)-volume densities si, S2 — si, • . . , Sd-k — Sd-k-i, respectively. Thus, 
by iterated application of intersection formulae for such processes (see \13\ Theorem 4.5.3]) 
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we find that Q'^]!'"'^''' '"^^ is proportional to the product si(s2 — si)(s3 — S2) • • • (sd-fc 



Sd-k-i) with some proportionahty constant ci only depending on the dimension para- 
meters 
obtain 



meters k and d as well as on the hyperplane measure A. Putting now C2 := ci/^j^^]., we 



d-k 

^{su...,sa.k,t) U^{Sj - Sj-i) 

P(/3i G (0, si), . . . , Pd-k G (sd-fc-i, Sd-k)) = — — = C2 — 

Sk,k 

with the convention that sq = 0. Differentiation implies that the joint density of the 
birth-times is C2t~^'^~^\ Since this integrates to 1, we must have C2 = {d — k)\, which 
completes the proof. □ 

To give a proof for general j £ {0, . . . ,k} we need the following fact. 

Lemma 6. Let d > 1, k £ {0, . . . , d — 1} , f : — >• M 6e non-negative and measurable, 
t > and < si < . . . < Sd-k < t- Then 

(8) E[/(Mp(;i)|/3i = su.. .,Pd-k = Sd-k] = E[/(F™^(^'^-'=))] . 

Proof. Since (Y{t))t>Q - regarded as a random process taking values in the measurable 
space of tessellations of ~ has the Markov property in the continuous time parameter 
t, it holds that 

E[/(MPg;^)|/3i = si,. . .,Pd-k = Sd-k] = ^[f{Mvf^^)\Pd-k = Sd-k]. 

Next, Lemma |4j tells us that the conditional distribution of MP^*]^, given its last birth 

time (id-k = Sd-fc, is the same as the distribution of F^^"'"*'**''"''^. Formally, this is exactly 
([s]), which completes the argument. □ 

3.2 Proof for general j 

The relationship between MP^*^. and MP^*q can be described by 

(9) E/(MPg) = [Ey,(MpWj]-^E[/(MP«)y,-(MpW)], 

where / is a non- negative measurable function on This is a direct consequence of 
Neveu's exchange formula [13, Theorem 3.4.5] and can be shown similarly to Equations 
(8) and (10) in [12j. In particular for j = k we have that 

(10) E/(MPi*i) = [Eyfc(MPl*J,)]-' E[/(MPgJ,) T4(Mp1*J,)]; 

see also |3]. Using ( [lO| ) with f{-)Vj{-)Vk{-)~^ instead of / there, we find that 

(11) E[/(MpW)y,(MpW)F,(MP«)-i] = [Ey,(MPi*i)]-^E[/(MP«)y,(MPg))]. 
Combining now ^ with ([n]), we see that MP^*^. and MPj,*], are related by 



:i2) E/(MPg) = 1^^^^^ E[/(MPW ) F,(MPgi) y.(MPgi)-^l 

EVj(MPj;.n) 



Let us fix < ri < . . . < rd-k < t and apply (12) with 



/ = G (0,ri), . . .,/3d-k G ird-k-i,rd-k)} 
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to obtain 

>(*) \ r- fn ^\ a. . f^AT^it) ■ 



P(/3i(MP^*;.) G (0,ri), . . .,(3d^k{MP^k!j) ^ ird-k-i,rd-k)) 



(13) = ^^'';^"^^;°; E[i{/3i(Mpgi) e (o,n),... 



...,/3,_fc(MPgi) G (r,_fc_i,rrf_,)}y,(MPl*i)yfc(MPgi)-i]. 
Conditioning on the birth-times and using Lemma [5] yield 
E[l{/3i(MpW) G (0,ri),...,/3rf_fc(MpW) G (rd_,_i, r,_fc)} y,(MPgi) yfc(MPgi)-i] 

= / / . . . / E[y,(MPl*i) yfc(Mp(*i)-i|/3i = si,..., Hd-k = Sd-k] 

Jo Jri Jra-k-l 

{d-k)\ ^ 
^ fd-k ^^d-k ■ ■ ■ ds2dsi. 

Lemma [ojimpUes that the joint conditional distribution of (Vfc, Vj) of the V^-weighted typ- 
ical A;-dimensional maximal polytope MP^*^^ of y(t), given its birth-times (/3i, . . . , Pd-k) = 
(si, . . . , Sd-k), only depends on the last birth-time Pd-k = Sd-k and equals the joint dis- 
tribution of (Vfc, Vj) of the Vfc-weighted typical fc-face F^^"'"^*''"*'' in a Poisson hyperplane 
tessellation with intensity measure Sd-k^- Whence, due to the scaling property ^ and 
the homogeneity of the intrinsic volume Vj we infer that 

E[^,(MPgi) 14(MPl*i)-i|/3i =si,..., f3d-k = Sd-k] 



k,k 



PHT{l)._n 
k,k ) J 



=43EK(Frr")i'*(Frr")-'] 

k— j 

where C3 := E [yj(F^^"'"^^^) ^^(F^^'^^"^'')"-'^] is a constant only depending on j, k, d and A. 
So, 

E[l{/3i(Mp(*h G (0,ri),...,/3d_fc(MPi*h G {Td-k-i.rd-k)}V,{MvfAVu{Mvf^r^] 



(14) 



Jri 



rd-k-i 



k-j {d-k)\ 
cssa-i ^a-k "^^d-k ■ ■ ■ dsadsi. 



Moreover, according to Lemmajsjb), there is another constant C4, only depending on the 
dimension parameters j, k, d and on the hyperplane measure A, such that 



(15) 



Eyfc(MPgi,) 
E^i(MPl*;,) 



C4 in- 



putting C5 := C3C4 and combining (13) with (14) and (15), we arrive at 



P(/3i(Mp(*)) G (0,ri),...,/3rf_fc(MpW) G (r,_fc_i, r,_fc)) 



Jri 



rd-k 



C5 (d - A;) ! 



''"d-fc-l 



k—j 
^d-k 



dsd-k ■ ■ .ds2dsi. 
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This shows at first that the birth-times of MP^*^- have a joint density, which is given by 

(si, . . . , Sd-^k) ^ c^{d- k)\ ^^1{0 < si < . . . < Sd-k < t}. 

Since this must integrate to one, we must have C5 = (d — j)/{d — k). This finally proves 
that 

k—j 

(si, . . . , Sd-k) ^{d-j){d-k- 1)! ^1{0 < si < . . . < Sd-k < t} 
is the joint birth-time density of MP^*^.. □ 

4 Proof of Theorem [2] 

We show the formula only for the case of the length- weighted typical maximal segment, 
the argument for typical maximal segment is similar. The proof follows the idea of the 



corresponding special cases d = 2 and d = 3 in (TJHTIIIH]. The key will be Equation (19) 
below, which is simple in the planar case and which has been established in [18] for d = 3 
in a long and intricate proof using other methods (parts of which seem to be restricted 
to small space dimensions). Here, we will employ the power of the distributional identity 
in Lemma H] combined with the result of Theorem [H 

We recall that the length-weighted typical maximal segment of Y{t) is the intersec- 
tion of d — 1 maximal polytopes of dimension d — 1, whose birth-times are denoted by 
/3i, . . . ,/3rf_i. Furthermore, let I and be the length and the direction of the segment, 
respectively, where (/9 G is the unit vector parallel to the segment. The distribution 
of Lf is given by a probability measure R on S'^^; cf. [4, Theorem 1] or [13, Theorem 
4.4.8]. However, the exact form of R is irrelevant for our purposes. 

First, by conditioning, we find that 

Pi,i(n) = / / ... / P{N = n\£ = x,ip = u,/3i = si,. . . ,/3d-i = Sd-i) 

Jo Jo Jo Js'l^'^ Jo 

xPe,l3i,...,Pa_^\v=u{x, si,..., Sd-i) dxR(du) dsi . . . ds^.i, 

where is the number of internal vertices and Pe,(3i,...,f3^_i\ip=u is the conditional joint 
density of length i and birth-times . . . ,/3d-i, given the direction u of the segment. 
Two comments are in order. At first, we notice that the joint density p^,(^,/3i,...,/3^_j. of 
length, direction and the d — k birth-times does not need to exist since R is not assumed 
to be absolutely continuous with respect to some reference measure on This causes 

the extra integration with respect to R and forces us to condition on (/?. It is also worth 
pointing out that the order of integration in the above expression cannot be changed 
because of the dependencies among the involved random variables. 

We now determine the exact expression for P£,/3i,...,/3^_i|(p=M- To start with, write 

^gs| P£,/3l,...A-ll¥'=«(^' si,---, Sd-l) 

= P£|</p=«,/3i=si,-A-i=Sd-i(^)^'/3i,-A-ilv=«(«l'- ■ ■ 'S<i-l), 

where Pe\^=u,[ii=si,...,fSa-i=Sd-i the conditional density of £ given that ip = u, (3i = 
si, . . . , Pd-i = Sd-i and where J'/3i,...,/3^_;^|(p=u is the conditional joint density of the birth 
times f3i, . . . , Pd-i given the direction if = u of the segment. Now, Lemma [6] applied with 
k = 1 shows that 

Pe\v^=u,^i=si,...,0a-i=sd-i (^) = P£\v:>=u,pd-i=Sd-i (^) 



9 



and that this is the same as the length density of the length-weighted typical edge of the 
Poisson hyperplane tessellation PHT(srf_i). It is well known (see Example 1.4 in pQ) that 
this length is gamma (Erlang) distributed with parameter (2, A{{u))sd~i), whence 

(17) P%=.A-i,... A_.=s,_, (x) = A((n))2.iixe-^«"))^-^^ 

Moreover, Lemma [6] also implies that the birth-times are independent of the direction of 
the segment, i.e., 

(18) p^,,...,^^_^|^=„(si,...,Srf_i) =p/3,,...,/3^_^(si,...,Srf_i) = -^^^rfTT^ 

(in fact the direction of F^^"'"^''^ does not depend on the intensity s; see fi|). Inserting 



(17) and (18) into (16), we arrive at 

(19) P.A,...A-il^=«(^> • • • > ^d-i) = ^^^A((n))2.2_^^e-^(<"»^-^^ 

It remains to determine the conditional probability 

p(Ar = n\£ = x,ip = u,/3i = si,.. . = Srf_i). 

For this, we first notice that for space dimensions d > 3, the length-weighted typical 
maximal segment can already have internal vertices at the time of its birth. Sup- 
pose now that the length i = x, the direction ip = u and corresponding birth-times 
/3i = si, . . . ,/3d-i = Sd-i with < si < . . . < Sd-i < t are given. We now gradually 
reconstruct the internal structure of the segment and exploit the intersection property 
of STIT tessellations. From time si until time ,52, behind the first {d — l)-dimensional 
maximal polytope. Mi say, a number of internal vertices appear, which are Poisson dis- 
tributed with parameter A{{u))x{s2 — si) since the region behind Mi can be further 
subdivided during the construction. Next, from time S2 until time S3, behind the first 
(d — l)-dimensional maximal polytope Mi there is another such polytope M2 so that 
internal vertices can be created by further subdivision within M^ n M^ and M^j'' n M^, 
where M.f stands for the half-space behind Mi (i.e., the last {d— l)-dimensional maximal 
polytope Md-i appears in M^). Thus, we obtain in the time interval (52,53) a Poisson 
distributed number of internal vertices with parameter 2A((n))a;(s3 — S2). In general, 
from time Sk until time Sk+i, 1 < k < d — 2, behind Mi we have k — 1 further (d — 1)- 
dimensional maximal polytopes M2 , . . . , M^ diving Mf into k regions so that by further 
subdivision a Poisson distributed number of internal vertices appears, whose parameter is 
kA{{u))x{sk+i — Sk)- During the last time interval {sd-i,t), we have not only to consider 
the number of internal vertices appearing within the d — 2 regions behind Mi, but also 
the number of internal vertices which arises in the two cells in M^ on the left and on the 
right of the segment (i.e., in the two regions separated by Md-i). This leads to a Pois- 
son distribution with parameter d ■ A({u))x{t — Sd-i)- Adding all the above independent 
Poisson distributed numbers, we obtain a Poisson distributed number of internal vertices 
with parameter A{{u))x{d ■ t — 2sd-i — Sd-2 — ... — si). Formally, this means that 

P(iV = = x,99 = u, /3i = si, . . . ,^d_i = Sd-i) 

_ {M{u))x{d • t - 2Sd-l - Sd-2 ^l))" -A(M)x.(d.f-2 



Thus, 



Jo Jo Jsl'^ Jo 



Pi,i(n)= / r'- ri_l ^^A((-))^«Li^e-'^(<"^^^-^^ 



d-i ./n t 



X 



{A{{u))x{d ■ t - 2sd-i - Sd-2 - ... - si))"" 
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Integrating now with respect to x, we see that all terms A{{u)) cancel out and we arrive 
at the desired formula for pi^i(n). 

The value of pi.i (n) is necessarily independent of t since the number of internal vertices 
does not change when the tessellation is rescaled. The latter, however, is because of Q 
the same as a time change, so that pi^i(n) must be independent of t. 

The formula for the mean number of internal vertices of the length-weighted typical 

oo 

maximal segment can be obtained from ^^npi^i(n) by a straight forward integration 

n=0 

procedure. □ 
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